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Abstract —  This  paper  reviews  recent  advances  in  designing 
unimodular  sequences  with  good  auto/cross  correlation 
properties  along  with  a  new  approach  that  emphasizes  on 
independent  receiver  design.  The  general  prohlem  is  to  design 
single  or  multiple  sequences  with  constant  modulus  in  the  time- 
domain  such  that  their  respective  matched  Alter  outputs  ideally 
resemble  delta  functions  and  the  cross-matched  filter  outputs  are 
zeros.  In  this  context  CAN  (cyclic-algorithm  new)  and  WeCAN 
(Weighted  CAN)  have  been  proposed  for  designing  such 
sequences  with  good  auto-cross  correlation  properties.  In  this 
paper,  the  equivalence  of  the  CAN  algorithms  and  the  classic 
Gerchberg-Saxton  (GS)  algorithm  involving  the  sequential 
magnitude  substitution  operations  in  the  time  and  frequency 
domain  is  demonstrated. 

The  design  of  unimodular  sequences  is  further  generalized 
here  by  considering  the  receiver  design  to  be  more  general  than 
the  respective  matched  input  sequences.  The  receiver  design  is 
carried  out  by  taking  care  of  the  desired  output  requirements 
and  the  freedom  present  at  the  input  can  be  used  to  further 
minimize  the  output  side-lobe  level. 

I.  Introduction 

In  a  series  of  excellent  papers,  Stoica  et.  al.  has  recently 
addressed  the  problem  of  designing  unimodular  sequences  with 
good  auto/cross  correlation  properties  [1,2,  3].  The  problem  is 
to  design  single  or  multiple  sequences  with  constant  modulus 
in  the  time-domain  such  that  their  respective  matched  fdter 
outputs  ideally  resemble  delta  functions  while  the  cross- 
matched  fdter  outputs  are  zeros.  Since  the  matched  fdter 
outputs  represent  auto  and  cross-correlations  of  the  original 
sequences,  the  problem  is  to  design  constant  modulus  time- 
series  whose  auto  correlations  approximate  the  delta  function, 
and  the  cross  correlations  approximate  the  all-zero  function.  In 
this  context  CAN  (cyclic -algorithm  new)  and  WeCAN 
(Weighted  CAN)  have  been  proposed  for  designing  such 
sequences  with  good  auto-cross  correlation  properties  [1,  3]. 
In  this  correspondence,  the  equivalence  of  the  CAN  algorithm 
and  the  classic  Gerchberg-Saxton  algorithm  is  demonstrated 
[4].  In  the  multichannel  case,  the  CANAVeCAN  algorithms 
coincide  with  the  Gerchberg-Saxton  algorithm  in  the  time 
domain,  and  they  generalize  into  the  frequency  domain  in  an 
interesting  manner. 

The  design  of  unimodular  sequences  is  further  generalized 
here  by  considering  the  receiver  design  to  be  something  other 


than  their  respective  matched  fdters.  The  receivers  need  not 
have  the  unimodular  property,  and  instead  they  can  emphasize 
on  the  desired  impulse-like  properties.  This  approach  gives 
rise  to  more  efficient  receiver  design  that  realize  the  desired 
output  sequences  and  the  freedom  present  at  the  input  can  be 
used  to  accommodate  additional  constraints. 

II.  Gerchberg-Saxton  Algorithm 

The  Gerchberg-Saxton  algorithm  considers  the  problem  of 
reconstructing  a  time  function  with  known  partial  magnitude 
information  both  in  the  time  and  frequency  domain.  Thus  let 

/(f)  =  a(f)e^'®"'  -fA F{co)  =  M(ffl)/'"<'“>  (1) 

represent  a  Fourier  transform  pair.  Given  the  information 

a{t)  >0,  0<t  <T,  and  M(®)  >0,  |(a|  <  ,  (2) 

both  in  time  and  frequency  domains,  the  problem  is  to 
reconstruct  the  entire  / (f)  in  some  optimal  fashion.  It  is  well 
known  that  the  set  of  all  signals  that  have  a  prescribed  value  in 
the  interval  (0,  T)  forms  a  closed  convex  set  C .  Convex  sets 
have  the  property  that  for  any  point  /(f)  outside  C,  there 
exists  a  unique  nearest  neighbor  Pf{t)  in  C  such  that  [5]. 

||/(f)-P/(f)||<||/(f)-g(f)||,  forallg(f)eC.  (3) 

Although  the  set  of  all  signals  with  the  given 
magnitude  transform  M{co)  over  a  prescribed  bandwidth  Q 
do  not  form  a  closed  convex  set  [5],  nevertheless,  it  is  possible 
to  assign  to  every  arbitrary  signal  / (f)  that  is  outside  this  set  a 

“nearest  neighbor”  signal  P^f{t)  that  belongs  to  such 

that  there  exists  no  other  signal  g(f)  e  for  which 

||/(0-g(0||<||/(0-/«/(0||  (4) 

is  satisfied.  In  the  case  of  closed  convex  sets  the  above 
operator  is  the  projection  operator  and  the  “nearest  neighbor” 
is  unique  [5].  In  the  case  of  the  above  magnitude  substitution 
operator,  the  uniqueness  property  is  not  preserved  for  the 
“nearest  neighbor”  and  in  general  only  the  inequality  constraint 

||/(0-/«/(0||^||/(0-g(0||,  forallg(f)eC„  (5) 
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is  satisfied.  Given  an  arbitrary /(t) ,  to  determine  P^f{t) ,  we 
need  to  generate  g{t)&C^  such  that  ||/(t)-g(t)f  is 
minimum  over  Q  .  With  / {t)  <->  F{(a)  = 


1  F{o}), 


ft)  e  Q 
ft)eQ' 


gives  the  magnitude  substitution  operator  [4,  6,  7]. 


(6) 


Similarly,  any  prescribed  temporal  magnitude  function 
a{t)  >0  as  in  (2),  or  in  particular,  a  constant  envelope  signal 

/ (t)  can  be  expressed  as  / (t)  =  where  ^  is  a  suitable 

positive  constant  that  can  be  used  to  maintain  a  prescribed 
energy  level  for  f  {f) .  It  is  easy  to  see  that  signals  with  the 

same  constant  envelope  do  not  form  a  convex  set;  however 
their  behavior  is  similar  to  the  signal  set  with  given  Fourier 
transfonn  magnitude.  Flence  if  denotes  the  set  of  functions 

{g(f)}  that  maintain  a  constant  envelope  level  A,  although 
is  not  a  convex  set,  it  is  possible  to  define  an  operator 
that  assigns  to  every  arbitrary  function  f  {t)  a  nearest  neighbor 
belongs  to  as  in  (5).  Thus  with 

/ (t)  =  b(t)  ,  we  have 


^./(O 


{AeF^‘\  te(0, 7) 

[  / (t),  otherwise 


(7) 


where  the  interval  (0,  T)  represents  the  duration  over  which 
the  constant  envelope  property  is  to  be  maintained.  The  error 
between  f^{t)  and  its  constant  envelope  version  P^fi^{t)  or 

the  transform  magnitude  substituted  version  of  (t)  reduces 
as  A:  — >  00  [6]. 


More  interestingly,  the  iteration  that  combines  both  the 
Fourier  transform  substitution  operation  as  well  as  the  constant 
envelope  substitution  operation  has  the  desired  error  reduction 
properties.  Thus  with 


L.At)=PAPMm=p,g,{t)  (8) 

where  we  define  g^it)  =  P/^ftif)  ■  ^^e  error 

reduction  property  [4] 

=  ||/a-+i  -  II  =  ||/tu  -  PMfk.i  II  ^  llA  -  -Pm A II  =  A  (9) 


is  maintained  since  (see  Fig.  1) 

A+i  ~  iia+i  ~7^a+i|| — iia+i  ~  A)  a| 

=  |AAA-7’mA|N||A-A.A||  =  A 


(10) 


In  (8),  /j,(t)  represents  the  constant  envelope  signal, 
whereas  gj(I)  preserves  the  desired  Fourier  transform 

magnitude  function.  Furthermore,  the  above  error  reduction 
property  is  maintained  also  if  the  above  iteration  includes 
convex  projection  operators  P  such  as  in  the  iteration 


f,At)  =  pp,p,Mt).  (11) 


Here,  for  example,  P  can  represent  the  projection  operator  that 
tmncates  the  constant  modular  signal  P^P^fi^{t)  to  the  interval 

(0,  T)  obtained  by  substituting  it  to  zero  outside  that  interval. 
The  iterative  algorithm  (11)  also  satisfies  the  error  reduction 
property  as  in  (9)  -  (10)  [7].  Equations  (8)  -  (11)  represent  the 
classic  Gerchberg-Saxton  algorithm  [4]. 


Fig.  1  Error  reduction  property  for  magnitude  substitution  operators. 

Observe  that  the  iteration  in  (8)  is  an  excellent  way  to 
prescribe  a  given  Fourier  transform  magnitude  while 
maintaining  a  constant  envelope  in  the  time  domain,  and 
together  they  result  in  iterations  that  “get  closer”  in  the  sense 
of  (9).  From  (8),  /j,^j(t)  satisfies  the  constant  modulus 
property  while  only  approximating  the  given  magnitude 
transform  M((a) ;  whereas  g*  (t)  =  A, A+i  (0  possesses  the 
given  magnitude  transform  while  only  approximating  the 
constant  modulus  property.  Nevertheless,  after  a  large  number 
of  iterations,  as  (9)  shows,  the  difference  between  these  signals 
gets  smaller.  Notice  that  since  (9)  -  (10)  do  not  guarantee 
convergence  properties  for  /j^(t),  and  in  general  there  is  no 
unique  solution.  The  final  solution  is  sensitive  to  the  starting 
point  f^(t) ,  and  a  careful  selection  of  the  initial  solution  based 
on  other  factors  is  essential  [4,  7,  8,  9]. 

In  particular,  the  discrete  version  of  the  Gerchberg-Saxton 
algorithm  described  above  can  be  summarized  as  follows  [10]: 


Let  [flY^  ™  A'-point  DFT  pair  at 

CO,.  =  iTtk/N ,  k=\^N ,  whose  magnitude  values  \f.\  =  a, 
and  17^^!  =  ^^,  are  specified  both  in  the  time  and  frequency 
domain.  Here  a,  =  a  for  example  will  refer  to  the  unimodular 
case.  Perform  the  A'-point  DFT  of  the  given  sequence  and 
substitute  the  frequency  magnitude  sequence  to  ,  then 
perform  the  inverse  DFT  and  substitute  the  temporal 
magnitude  values  to  a,  and  repeat  the  procedure.  For 
example,  Wikipedia  summarizes  the  above  algorithm  as  [10]: 


Gerchberg-Saxton  Algorithm ( Source ,  Target , 
Retrieved_Phase ) 

A  =  IFT (Target) 

while  error  criterion  is  not  satified 
B  =  Amplitude ( Source )  *  exp (i*Phase (A) ) 
C  =  FT (B) 

D  =  Amplitude (Target)  *  exp (i*Phase (C) ) 
A  =  IFT (D) 

end  while 

Retrieved_Phase  =  Phase (A) 
end  Gerchberg-Saxton  Algorithm 


(12) 


The  error  reduction  property  in  (9)  is  guaranteed  in  this 
case  also.  Interestingly  for  the  7/-point  discrete  case,  the  above 
algorithm  implicitly  contains  the  temporal  limiting  projection 
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operator  P  described  in  (11)  and  hence  it  represents  the 
iteration  in  (11). 

III.  Unimodular  Sequences  uisng  CAN  And  WeCAN 


optimization  problem  A  =  min||5(c(;)-/||  is  more 

meaningful.  It  also  involves  the  difficult  task  of  minimizing 
the  error  over  all  frequencies.  However,  the  relaxed  criterion 


In  the  single  channel  case,  as  Fig.  2  shows,  the  problem  is 
to  design  a  sequence  such  that  |xj  =  1  for  ;  =  1->A, 

and  their  matched  filter  outputs  approximate  an  ideal  delta 
function  for  pulse  compression  and  sidelobe  suppression. 
Since  the  matched  filter  outputs  Ir,  1™'  are  the 

autocorrelations  of  the  input  sequence,  we  have 

h  =  Z  =^-k’  ^  =  0>  ±1>  ±(A^-1)  ■ 

/=A'+1 


Fig.  2  Matched  filter. 

Hence  for  pulse  compression  and  side  lobe  suppression, 
ideally  we  must  have  a  delta  function  at  the  output,  i.e.,  =1 , 

=  0,  /t  >  1  or  in  this  context,  minimization  of  the  total 
sidelobe  error  [3] 

e  =  {r^-\f+2^\r,\^  (14) 

k=l 

has  been  suggested  as  a  suitable  criterion.  The  minimization  of 
the  above  error  is  closely  related  to  the  minimization  of  the 
integrated  sidelobe  level  [ISL) 

ISL  =  Y\r,\\  (15) 

k=\ 

or  the  maximization  of  the  merit  factor  {MP^  MF  =  jllSL 

proposed  in  the  literature  [11,12].  It  is  easy  to  relate  the  error 
in  (14)  or  the  ISL  in  (15)  to  the  classic  discrete-time  Fourier 

...  N 

transfonn  of  the  sequence  |x,-|  given  by  X(rE»)  = 

i=i 

Using  (13),  the  corresponding  spectmm  S((a)  is  given  by 

S(o))  =  \X{co)f  =  ^  r,  (16) 


A. 


1 

2N 


IN 

ZI|A(®„)-/| 


<y.. 


2  Tin 
2N 


(18) 


obtained  by  averaging  the  spectral  error  only  at  2N  equally 
spaced  discrete  points  is  more  appealing  since  it  is  also  equal 
to  the  total  sidelobe  error  s  in  (14).  Using  (16),  we  also  have 


and  clearly  (19)  is  minimized  by  letting  [1] 

X(cu„)  =  e"'"",  n  =  \^2N .  (20) 

X((e)^J  represent  the  DFT  coefficients  at  frequencies 
=2Knl2N,  and  constraint  (20)  is  implementation  of  the 
unity  magnitude  substitution  requirement  in  the  discrete 
frequency  domain  as  in  (12)  .  Using  the  2N  y-2N  DFT  matrix 

D  whose  (i,  j)"'  element  is  given  by  £)  = _ I _ ^jinikUN  ^ 

we  can  rewrite  the  2N  equations  in  (20)  compactly  as 

D*A  =  ,  •  ■  ■  J  =  V,  or  X=DV  (21) 

where  the  2N  x  1  vector  X 

A  =  [x,,  Xj,  ■■■Xjv,  0  ■■■  O]’^,  |x,.|  =  l  (22) 


represents  the  extended  data  vector  with  x,. ,  /  =  1  — >  A 
representing  the  unimodular  entries.  Equations  (21)  -  (22) 
suggest  that  for  the  above  minimization  problem,  the  2N  - 
point  DFT  of  the  constant  modulus  data  vector  X  in  (22)  in 
the  temporal  domain  should  generate  another  constant  modulus 
vector  V  in  the  frequency  domain  and  the  CANAVeCAN 
algorithms  perform  these  operations  by  employing  the 
magnitude  substitution  operation  both  in  the  frequency  domain 
and  the  time  domain  sequentially  till  the  desired  accuracy  is 
achieved  [1,  2,  3].  But  this  is  the  same  as  a  Gerchberg- Saxton 
algorithm  described  in  (8)  -  (12).  The  CAN  algorithm  in  (21)  - 
(22)  resets  the  second  half  of  the  temporal  coefficients  to  zero 
in  (22)  at  every  stage,  and  it  corresponds  to  the  projection 
operator  P  associated  with  the  time-limiting  operation  in  (1 1). 


and  hence  when  the  ideal  conditions  where  is  delta  function 
are  substituted  into  (16)  we  obtain  the  new  requirement 

(17) 


IV.  Multi-Channel  Case 

The  multichannel  version  of  the  CAN  approach  also  uses 
the  Gerchberg-Saxton  algorithm  in  the  time  domain  while 
generalizing  it  in  frequency  domain.  The  NxM  matrix  [2,  3] 


for  all  0)  in  the  frequency  domain.  Eq.  (17)  represents  the 
well  known  result  that  to  achieve  excellent  pulse  compression, 
the  input  sequence  must  exhibit  flat  spectrum.  Since  the 
requirement  in  (17)  is  impossible  to  maintain  over  all 
frequencies  by  nontrivial  finite  sequences,  to  generate 
unimodular  finite  length  sequences  that  satisfy  (17),  the 


X  =  [xi,  x^,  ■■■  x„]  = 


T 

T2 


T 

Yn 


(23) 
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containing  the  M  temporal  sequences 

i  =  |x,.^.|  =  l  (24) 

represent  M  orthonormal  sequences  each  with  N  unimodular 
samples,  where  the  «'*  row 

/„  =  \.x^{n),  x^{n),  ■■■  x^{n)\  (25) 

in  (23)  corresponds  to  the  M  channel  outputs  at  the  «'*  time 
instant.  Under  stationary  assumptions,  cross-correlations  of 
different  channel  waveforms  give  rise  to  the  auto/cross¬ 
correlation  coefficients 

^i.  j  (^)  =  S  (« -  ^)  • 

;i=l 

Using  (25),  the  MxM  auto/cross-correlation  coefficient 

matrices  generated  using  (26)  can  be  compactly  expressed  also 
using  the  M  x  1  column  vectors  as 

=  ^  =  0^A^-1  (27) 

n=\ 

where  represents  the  MxM  cross-covariance  matrix 
corresponding  to  lag  k .  For  orthonormal  sequences  we  require 

R„=I,  R,=0,  ^>1.  (28) 

Following  (14)-(15),  if  we  define  the  total  sidelobe  error 
criterion 

^  =  ||R„-lf+2X||R,f  ,  (29) 

*■=1 


in  the  frequency  domain.  Once  again,  Eq.  (33)  can  be 
compactly  expressed  as 


'  z^(®.)  ■ 

T 

a, 

Z"(®,) 

=  D 

X 

ii 

_Z^(®2A-). 

T 

2JV_ 

where  X  is  defined  as  in  (23)  with  \x.  ;|  =  1  ■  Observe  that  the 

normalization  of  a,,  in  the  frequency  domain  in  (33)  -  (34)  is  a 

generalization  of  the  scalar  version  of  the  Gerchberg- Saxton 
algorithm  since 

IKII  =ZK(«r)|  =M{ca,)  (35) 

m=\ 

represents  the  sum  of  the  magnitudes  of  the  various  channel 
components  of  the  DFT  coefficients.  Following  (33)  -  (34), 
the  multichannel  version  of  the  CAN  algorithm  performs  the 
magnitude  substitution  operation  sequentially  both  in  the 
frequency  domain  ( ||^  =  1 )  and  in  the  time  domain 

( |x.  .  I  =  1 )  as  in  the  Gerchberg-Saxton  algorithm  in  (12). 

Fig.  3  shows  the  matched  filter  receiver  outputs  of  a  ten 
channel  unimodular  sequence  of  length  200  generated  using 
the  CAN  algorithm  in  (33)  -  (34).  Both  Fig.  3(a)  and  (b) 
contain  ten  auto-correlation  sequences  and  90  cross-correlation 
sequences.  Notice  that  the  auto/cross  correlation  terms  are  at 
about  20  dB  below  the  peak  term.  The  sidelobe  levels  here  can 
be  further  improved  by  considering  various  weighting 
sequences  as  well  [2,  3]. 


N 

then,  as  before  with  =  we  have  the  power 

1=1  ~ 

spectral  density  matrix 

S(cu)  =  Z(£u)Z*(cu)=  Z  Rre-'*'"  (30) 

yt=-(V-l) 

SO  that  once  again  the  desired  error  criterion  in  (28)  can  be 
expressed  as  S((y)  =  I .  As  in  (18),  the  relaxed  criterion 


Delayr 


1  IN 


03,.  =- 


iTik 

2N 


(31) 


is  more  useful  in  this  context  since  it  relates  the  point  power 
spectral  error  with  the  total  sidelobe  error  in  (29).  Expanding 
(3 1)  using  (30)  we  also  obtain  [3] 


SO  that  minimization  of  s  is  achieved  by  setting 

Z.{03^')  =  a„  ||a,f  =1,  A:  =  1^2A  (33) 


(a)  Front  view 


(b)  Side  view 


in  (32).  Notice  that  a^.  in  (33)  are  the  2N-point  DFT  vectors, 
and  their  normalization  requirement  is  a  magnitude  constraint 


Fig.  3  A  ten  channel  unimodular  sequence  design  using  the 
multichannel  CAN  algorithm  in  (33)  -  (34).  Each  figure  contains 
ten  autocorrelation  sequences  and  90  cross-correlation  sequences. 
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V.  General  Unimodular  Sequences 

The  design  of  unimodular  sequences  can  be  further 
generalized  by  considering  the  receiver  design  to  be 
independent  from  the  actual  matched  fdter  sequences.  This 
approach  is  useful  in  many  applications  including  radar,  where 
the  receiver  does  not  have  to  be  unimodular  and  can  be 
something  other  than  the  matched  fdter  sequence.  In  this 
context,  the  problem  is  to  design  M  input  unimodular 
sequences  x.,  i  =  \^M,  |x;^.|  =  l  each  of  length  N  as  in 

(23),  and  their  corresponding  receivers 

=[\o  ^2’  /  =  1^M  (36) 


such  that  the  receiver  outputs  should  satisfy  the  “desired” 
output  sequence  properties.  Thus  with 

^i,M)  =  '^xXk)KXn-k  +  \),  n  =  \^N  +  K-\  (37) 

representing  the  output  at  the  m'*  receiver  due  to  the  i"' 
sequence  at  time  instant  n  ,  then  ideally,  we  require  (Fig.  4) 


i  =  m 


0, 


(38) 


i^m 


where  corresponds  to  the  desired  peak  term  at  the  output. 

In  general,  the  receiver  length  AT  is  a  free  parameter  left  to 
other  design  considerations. 


Fig.  4  Generalized  unimodular  sequences. 

One  approach  in  this  context  is  to  use  a  criterion  similar  to 
the  merit  factor  (MF).  Thus  if  we  define 


'7,=- 


k, ■,■(«„)!' 


i=l^M 


(39) 


j^i  n=\ 


to  represent  the  ratio  of  the  receiver  output  peak  value  power  to 
the  total  undesired  sidelobe  power  at  each  receiver  output,  then 
the  receiver  design  problem  is  to  maximize  rj.,  ;  =  1  — >  M  . 

Observe  that  the  criterion  in  (39)  is  well  known  and  it  is  the 
same  as  the  approach  behind  the  classical  matched  filter.  Eq. 
(39)  can  be  expressed  in  matrix  fonn  as 


T].  =  '  ^  ,  1  =  1— >M  (40) 

h]A,h, 

where 

a,.=f,:oF,„+xf;f;>o  (41) 

represents  a  KxK  positive-definite  matrix,  and  F,.  represents 


the  {N  +  K-\)y.K 

matrix 

^  A(l) 

0 

0 

0 

x,.(2) 

x,(l) 

0 

x,-(2) 

0 

F,.  = 

x,(A) 

x,.(l) 

0 

0 

x,(A) 

x,.(l) 

0 

0 

x.(A-l) 

x,(A) 

x,.(A-l) 

.  0 

0 

0 

0 

x,(A) 

generated  from  the  i'*  input  sequence  x-  ■  In  (41),  F,.  ^ 

represents  the  {N  +  K  -  2)x  K  matrix  generated  from  F  with 
its  row  deleted.  Finally  in  (40),  represents  the  row 
vector  correspond  to  the  row  of  F,. .  Thus  in  general 

b„  =[0,  •••  0,  x,(iV),  x,(iV-l),  •••  x,(l),  0,  •••  0]  (43) 

for  n^>  N  ■  The  number  of  zeros  in  (43)  depends  on  the  value 
of  with  respect  to  N .  For  example,  with  n^=N  and 
Ai  =  ,  the  vector  is  the  same  as  the  flipped  version  of  x . 

in  (25).  Schwarz  inequality  applied  to  (40)  gives  the  optimum 
receiver  filter  to  be 


^,=A-‘£^.  (44) 

Observe  that  the  receivers  h.  in  (44)  are  obtained  non- 

iteratively  and  the  freedom  present  at  the  input  makes  the 
unimodular  constraint  trivial  to  implement  at  the  input. 
Although  each  receiver  design  depends  on  all  input  sequences 
through  the  matrix  A  in  (41),  nevertheless  they  can  be 

determined  in  an  uncoupled  manner  as  in  (44).  Furthermore, 
the  explicit  design  for  the  receiver  allows  freedom  in  deciding 
the  output  instant  at  which  the  peak  should  be  observed. 

This  freedom  allows  the  various  peak  outputs  to  be  time 
aligned,  a  desired  property  in  some  radar  applications. 

Fig.  5  shows  the  receiver  output  using  the  new  design 
approach  in  (44)  for  a  ten  channel  randomly  generated 
unimodular  sequences  of  length  200.  For  comparison 
purposes,  here  the  receiver  length  is  set  equal  to  200.  Notice 
that  the  auto/cross  sidelobe  terms  are  about  20  dB  below  the 
peak  value,  similar  to  that  in  Fig.  3.  This  represents 
performance  similar  to  the  CAN  sequences  obtained  using  (33) 
-  (34).  In  this  new  approach,  additional  optimization  can  be 
carried  out  at  the  input  to  further  reduce  the  sidelobe  level  and 
for  compression  in  the  Doppler  domain. 
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(b)  Side  view 


Fig.  5  A  ten  channel  unimodular  sequence  auto  and  cross  terms 
at  the  output  of  the  sidelobe  minimizer  receiver  in  (44).  (a)  Front 
view,  (b)  Side  view. 

VI.  Conclusions 

This  paper  reviews  recent  advances  in  designing 
unimodular  sequences  with  good  auto/cross  correlation 
properties  along  with  a  new  approach  that  emphasizes  on 
independent  receiver  design  by  minimizing  the  sidelobes. 
Since  the  matched  fdter  outputs  represent  auto  and  cross¬ 
correlations  of  the  original  sequences,  the  problem  is  to  design 
constant  modulus  time-series  whose  auto  correlations 
approximate  the  delta  function,  and  the  cross  correlations 
approximate  the  all-zero  function.  In  this  context,  CAN 
(cyclic-algorithm  new)  and  WeCAN  (Weighted  CAN)  have 
been  proposed  for  designing  such  sequences  with  good  auto¬ 
cross  correlation  properties.  The  equivalence  of  the 
CANAVeCAN  algorithms  and  the  classic  Gerchberg- Saxton 
(GS)  algorithm  involving  the  sequential  magnitude  substitution 
operations  in  the  time  and  frequency  domain  is  demonstrated 
here,  with  an  interesting  generalization  of  the  GS  algorithm  to 
the  multichannel  case. 

The  design  of  unimodular  sequences  is  further  generalized 
here  by  considering  the  receiver  design  to  be  independent  from 
the  actual  input  sequences.  The  receiver  design  takes  care  of 
the  sidelobe  requirements  and  the  freedom  present  at  the  input 
can  be  used  to  further  minimize  the  output  sidelobe  levels. 
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